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1. Introduction 



Restriction of a many body system to low dimension brings with it simplifications but 
also a cost - the need to pay much closer attention to the effects of both thermodynamic 
fluctuations and interparticle correlations. One-dimensional models have proved an 
extremely rich field for the exploration of subtle effects of electron-electron interactions. 
They have been found to have quite complex ground states brought about by the same 
reduced phase space which makes some of their properties soluble, and which causes 
finite temperature to play so dramatic a role. 

A colourful variety of quasi-one- and two-dimensional materials is now accessible to 
experiment, but this has not produced scientific consensus in all areas. Nine years of 
intensive research on high-temperature superconductors, the best known examples of 
quasi-two-dimensional systems, have failed to yield agreement on even the symmetry of 
their respective ground states. 

The study of quasi-one-dimensional conductors has been less intense, but they are 
interesting for similar reasons. All attempts to model these materials must make 
reference to the fact that they are not truly one-dimensional, but merely very, very 
anisotropic. The set of materials to which (TaSe^/ belongs have a common structure 
of an assembly of weakly coupled conducting chains. Their behaviour has generally 
been understood in terms of those properties of one (or three) dimensional systems 
which are held to be relevant to them. This interplay of dimensionality is not trivial; 
weak interchain coupling will act to stabilise states born of instabilities in the underlying 
one-dimensional structure, which would otherwise be destroyed at finite temperature by 
large thermodynamic fluctuations. In the limit of strong coupling both the instability 
and the fluctuations may be irrelevant. 

Many and diverse experimental techniques have been brought to bear on quasi- 
one-dimensional conductors in pursuit of insight into their properties and the low 
dimensional physics which these betray. The transition-metal tetrachalcogenide 
(TaSe^/ alone, in the thirteen years since its first synthesis [j], has been probed with 
neutron scattering |J, X-rays ||, low energy electrons ||, ultraviolet photons f|, and 
subjected to measurements of optical and electrical conductivity || [7|, |[ and magnetic 
susceptibility. |J. This material is interesting, since on cooling it displays a second 
order phase transition at the (relatively high) temperature of 263 K, from a highly 
conductive but extremely anisotropic "metallic" phase into a semi-conducting charge 



density wave ground state. While angle integrated photoemission [|1(], |TI| and inverse 
photoemision [[[2] experiments have been performed on such materials for some time, 
only very recently has angle resolved photoemission (ARPES) data become available. 
This offers the most direct means yet of understanding the changes which take place in 
the system's electronic properties. 
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The instability of one- dimensional metals against an insulating ground state was 
described by Peierls [13] in 1953 and the mean field theory [14] can be constructed 
in close analogy to the BCS theory of superconductivity. It is also well known that 
the fluctuations forced by the restricted phase space of a truly one-dimensional system 
prevent it from undergoing a phase transition at any finite temperature. The real 
systems, then, must exhibit an interplay of dimensional and thermodynamic effects, 
which manifest themselves in the strong and unusual temperature dependence of the 
"metallic" properties observed above the transition temperature. Early measurements 
of the magnetic susceptibility of (TaSe^/ H suggested that fluctuations of the charge 
density wave (CDW) order parameter were present for all temperatures above 263K, 
up to the limit of the compound's chemical stability (at about 430K). In order to 
understand the phase transition which takes place in this and other similar materials, 
it is necessary first to understand the role of fluctuations in the properties of a one- 
dimensional conductor. 

It has also been known for some time that the Landau Fermi liquid state is unstable 
against interparticle interaction in one dimension; the paradigm of the the Luttinger 
liquid has evolved to describe the properties of a number of abstract one- dimensional 
models (for an overview see Luttinger liquids have unusual correlation functions 

displaying separation of spin and charge degrees of freedom, and possess no stable 
single particle excitations at the Fermi energy ef- This implies that the customary 
discontinuity in zero temperature occupation number n(k) at the Fermi wavevector 
kp found in all Fermi liquids is absent, and that the electronic density of states will 
vanish near ep, with a power law behaviour determined by the strength of interparticle 
interactions. The decoupling of spin and charge degrees of freedom manifests itself in 
new structure in the spectral function of the system. 

A topic much discussed in recent years is whether some aspects of this well- 
established one-dimensional behavior can survive in higher dimensionalities. For 
example, it has been proposed that the normal state properties of high-temperature 
superconductors can be understood on the basis of this hypothesis [16], while others 
believe that any coupling between chains must destroy the Luttinger liquid. One 
approach to this problem is empirical: if evidence of Luttinger liquid behavior could 
be found for a real weakly-coupled-chain system this would constitute proof of the 
possibility of such behavior in higher dimensions. In this spirit, we ask whether 
(TaSe-^^I is, experimentally, a Luttinger liquid, as has sometimes been claimed. As 
a first step, we compare ARPES data with a more conventional interpretation, that of 
CDW fluctuations. In a later paper, we shall attempt a more systematic comparison 
of strongly correlated electron and CDW theories for the transport properties and the 
ARPES and core- level lineshapes of {TaSe^^I. 

In Sees. 2 through 4 of this paper we shall present the model used, its overall 
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properties, and the comparison of theory and experiment for ARPES lineshapes in 
(TaSe±)2l- The discussion is based on the premise that fluctuation effects born of the 
electron-phonon interaction dominate over those of correlation (interparticle interaction) 
for this system. It is found that such a treatment works well for (raSe?)^ . A critical 
comparison of this picture with other candidate descriptions is given in the conclusions 
of Sec. 5. 



2. The Model 



The problem of describing the combined effect of correlation and electron-phonon 
interaction in one dimension on an equal footing is an axiomatically hard one, 
since correlation in one dimension destroys the Fermi-Liquid picture on which the 
(perturbative) treatment of electron-phonon interaction is based, whilst the physics of 
the Peierls transition is dictated by 2k f 'backscattering' events, which cannot be treated 
within the usual scheme for electronic correlations. For a given real material, in a given 
temperature range, however, it may not be necessary to solve the general problem in 
order to understand the results of experiment. As there is no clear framework for relating 
the microscopic properties of materials to the parameters of any of the phenomenological 
models refered to in this paper, the best that can be done is to develop a model based 
on what is believed to be the relevant subset of physics for each material, and to test it 
against a variety of experiments. 

In order to provide a framework for our calculation we briefly review here some basic 
perturbative (mean field) results for electrons in a Peierls CDW System, and present 
a succinct derivation of a Green's Function for electrons suffering fluctuations of CDW 
order, in the spirit of the treatment of Lee, Rice and Anderson (hereafter referred to as 
LRA, 0). 

The starting point for any non-correlated theory of electron phonon interaction is 
the Frohlich Hamiltonian 

h = < k ) c Uk + Y.^iiHh + ~rf Y,9{q)ci +q c k u q (i) 

k q V L q,k 



u, = -7— (&J + &_,), 

where c\ and b q are (respectively) creation operators for electrons and phonons with 
dispersion e(k) and oo(q), u q is the Fourier transform of the lattice displacement, and 
g(q) the electron- ion coupling. L is the number of sites. For simplicity we will not credit 
the electrons with spin, and will describe the band by a free particle dispersion relation, 
linearised about kf. 
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Consideration of the linear response of this system |TJ| reveals that the phonon 
frequencies undergo drastic adjustment if the band is near half filling; the 2k f phonon is 
softened and its frequency vanishes at a well defined mean field transition temperature 
T c given by 

k B T c = —e~^v f kf, (2) 

7T 



(3) 



with 

A = JV„! 

h ^2k F 

where N Q is the density of states at the Fermi energy, and 7 ~ 0.5772 is Euler's constant. 
The overall picture of the important quantities in the mean field theory is presented in 
Fig. 1, where the gap A(T) and the phonon frequency u 2kF (T) are plotted on a phase 
diagram. 

The result in Eq. [| (the giant Kohn effect) may also be found in low order 
perturbation theory, in which case a length scale £0 emerges naturally from the 
calculation 



C(3) ~ 1.202 being the Riemann zeta function, and vf the Fermi velocity. Condensation 
of the softened 2kp phonon leads to a static lattice distortion < u 2kp >, and the opening 
of a gap A in the electronic spectrum. The interaction term in the Frohlich Hamiltonian 
(1) may then be approximated by a BCS-like mean field form 

H ( e f F) = £[ A *4-2A^ + Acl k+2kF c_ k ) (5) 



g(2 kp ) < u 2kF > 



Self consistent solutions of A and < u 2 k F > as a function of temperature may now be 
found for the Peierls insulator in the same way as for a BCS superconductor. Both 
systems are dominated by the same square root singularity in the electronic density of 
states at the edges of the gap. 

This picture of a metal-insulator transition in a purely one-dimensional metal is 
clearly not adequate as we expect fluctuations of the order parameter, when properly 
accounted for, to destroy the mean field solution at any finite temperature. Real systems, 
however, are not truly one-dimensional. {TaSe^) 2 I comprises parallel chains of tantalum 
atoms p0| , surrounded by approximately perpendicular rectangles of Se (Fig 2.). The 



iodine resides between chains. Overlapping d z orbitals on the tantalum chains form a 
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band along which conduction occurs. Whilst electronic transport across the chains is 
believed to be diffusive at the temperatures of interest here and thus does not give rise 
to coherent dispersion in the perpendicular directions, interchain interactions can act 
to stabilise the mean-field solution at some three-dimensional ordering temperature 
T C 3D considerably less than the mean field temperature T c . This then corresponds 
to the transition observed in real systems. The LRA model of a one-dimensional 
metal does not explicitly include interchain effects, but their relevance and the three 
dimensional ordering temperature emerge very naturally from their analysis. It is found 
that T 3D « T c /4. 

The physical content of Lee, Rice and Anderson's extension of the mean field picture 
is the realisation that the 2kf fluctuations of the softened lattice are slow on the timescale 
of electron dynamics and that they may therefore be approximated by a static disorder 
potential, the determination of which is then a separate (classical) problem. 

We start, then, from a natural generalisation of the BCS approximation to the 
electron phonon interaction 

He P = \^-Q C l>-Q C k> + ^Q c -k'+Q C -k'l (6) 

Q,k'>0 

where 

= -^=g{Q)u Q (7) 

are the components of the disorder potential and may be likened to an order parameter 
for fluctuations of some portion of the lattice. LRA then prescribe an equation of 
motion treatment of this hamiltonian, which after certain approximations generates the 
following relation for the electronic Green's Function 

G(k, k; iu} n )~ l = e{k) - iu n - ^Q^-q 77 n\ ~ " ( 8 ) 

The same relation may be written down immediately by using the analogy with static 
disorder simply to calculate the second order self energy correction for electrons scattered 
by the potential |TJ|. The Feynman diagram used is shown in Fig. 3. The box in 



the diagram is the phonon self-energy, or the charge-charge correlation function at 2k f. 
The deformed lattice approximation treats the ionic configuration as rigid (incapable of 
recoil), with the positions given by a thermal average. The absence of a frequency sum 
in the elecron self energy expresses the fact that the lattice and CDW fluctuations have 
been decoupled. 

We note in passing that it would in principle be possible to substitute a Luttinger 
Liquid Green's Function in the expression for the self energy, and so to treat the effects 
of lattice distortion on a correlated system to a similar level of approximation. To 
obtain and parameterise spectral functions from such a calculation would not, however, 
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be trivial. Similarly, refinements may be made to allow for lattice recoil. We shall not 
present these here, but will limit our discussion to the Green's function previously found 
by LRA. 

Evaluation of this expression then requires knowledge of the correlation function for 
the lattice fluctuations at a given temperature. A means of finding this was supplied 



by Scalapino, Sears and Ferrell [ I9| . They perform a functional minimalisation of a 
generalised Ginzburg-Landau free energy to obtain correlation functions for a fluctuating 
order parameter in a one-dimensional system, and obtain the following form 

< V(x)V{x') > = <^ 2 (T) >exp[-\x-x'\C 1 (T)]cos[2k F (x-x')]. (9) 

LRA substitute a free energy with parameters taken from the mean field (linear response) 
perturbative treatment of the ID Frohlich Hamiltonian : 



with 



F[* Q ] = a(T)\V Q \ 2 + 6(T)|vl/ Q | 4 + C (T)(Q - 2k F y\V Q \\ (10) 



T — T 



and 



6(r)=D [6o + (6d-6 1 )^: 



c(T)=D o e o (T), 



where Do is the (constant) density of states for the band, which is taken to have width 
2ep. We fix 6 and b\ to give the correct zero temperature value of the gap A in the 
electronic spectrum: b = and 



7C(3) (1.76) s 

01 — Oq- 



16tt 0.5 

The problem of determining < ^q(T) > and £(T) then reduces to that of finding 
the low lying energy levels of a particle moving in an anharmonic potential well, the 
shape of which is determined by the coefficients of the free energy 

H = - i -^ L ^-r + ^Hl*! 2 + Wl*| 4 (ii) 

This may solved numerically, or approximately using perturbation theory and 
asymptotic analysis. It is found that the coherence length £ increases steadily from 
its mean-field value at T c with reducing temperature, but increases very rapidly at 
a temperature approximately one quarter of T c . This implies that long-range order 
exists for temperatures below jT c , and interchain coupling stabilises the mean field 



solution. We will identify this temperature with the transition temperature T C 3D of a 
three-dimensional system and not attempt to treat interchain effects explicitly. The 
mean square value of the correlations increses approximately linearly with decreasing 
temperature, and takes on the role of a mean field gap below the three dimensional 
ordering temperature. 

The following parameterisation is accurate in the temperature range of experimental 
interest in the next section: 

AT 1 

r 1 (T)=^ 1 (T)(— --) (12) 

< *2 (T) >= _ ?_) _ l k = II 1 (13) 
K J b v T c 2 a' Jl-I- 

where a' = a(T)/T. 

We are now in a position to assemble the Green's Function for the system, with 
the (static) lattice fluctuations parameterised by the "gap" (squared) energy scale 
< ^ 2 (T) > and the "lifetime" energy scale £ 1 (T). 

G R (k, k-u)- l = uj- e(k) - [ dQ S{q )£Xl-v ( 14 ) 

J uj — e{kZQ) + to 

where S(Q) is a Lorentzian of width £ _1 centered on Q = 2kp, the sign in the 
denominator is chosen to give e(ktQ) ~ c(k), and < \I/ 2 > is found from the results of 



SSF 1 19). Evaluating the integral over Q, and dropping the second momentum index, 



we arrive at a result for the Green's function 



n ,, x_ g + e(fc) + iv F j 1 , . 

G R {k, u) - ^ _ <k)2 _ <r> +ivF ^ 1{u _ e{k)) ■ W 

This will form the basis for most of the subsequent analysis. 
3. Basic Properties of the Model 

Insight into the properties of the model outlined may be obtained by consideration of 
the imaginary part of the Green's function (15) derived above: 

A ^ u ) = TT2 77TX-2 Z , T ,2 ^12 , .. 2 ^-2r.. ^ia » ( 16 ) 



to 



2 - e(k) 2 - < ^ 2 >] 2 + vU~ 2 [u ~ <k)} 2 ' 



where the parameters < \E' 2 > and £ _1 have scale and temperature dependence given 
above. As observed by LRA, this may be integrated analytically to give an expression 
for the density of states. This is plotted for a system at temperature T = 300K= 
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with T c = 892K, T 3D = 263K and e F = 1.2eV in Fig. 4; the reason for this choice of 
parameters will be discussed in the light of photoemission data in a later section. 

The density of states shows clear evidence of a quasigap at all temperatures above 
T 3D observable for a system such as (TaSe^)2l ■ Spectral weight is still present at the 
Fermi energy (u = 0) at room temperature, but is greatly reduced. Traces of the square 
root singularity which will dominate the mean field solution are visible at the edges of 
the gap for temperatures approaching T C 3D . Sharp spikes and edges in the plot are a 
numerical artifact only. 

The LRA Green's function (15) above can clearly be seen to reduce to a BCS-type 
Green's function 

Gr ^ u) = u,*-e{k)*-<r>+*> ' (17) 

with < ^ 2 > taking on the role of a real gap in the limit where £ — > oo; this mean field 
solution is in fact the exact zero temperature limit of the LRA theory. At the level of 
approximation relevant to the experiments the model may be taken to posses a real gap 
below T 3D . We shall then proceed to describe the system by an LRA Green's function 
(15) above T 3D (263i^ for (TaSe^I), and by a BCS Greens function below T 3D . 

It is also possible to integrate the spectral function numerically over oo to obtain 
a result for n(k). (The integration my be performed analytically for the BCS-like 
expression below T 3D ). This is displayed in Fig. 5. The occupation number is clearly 
dominated by the presence of the quasigap, varying over a scale in k-space given by 
A/hvp- No Fermi step is present in the occupation number at an experimentally 
observable temperature, but the region of k-space over which n(k) undergoes most 
change becomes smaller with increasing temperature and decreasing size of the quasigap. 
Both the occupation number and the density of states illustrate the fact that the metallic 
(fluctuating CDW) phase of the model does not resemble a conventional Fermi liquid. 

Consideration of the denominator of the spectral function at hp shows that spectral 
weight is concentrated in one peak provided that the quantity < ip 2 > —v F C,~ 2 /4 is 
negative. This gap like parameter will change sign at a temperature of similar magnitude 
to, but in general different from, the mean field temperature T c . On cooling it becomes 
positive, and spectral weight is split into two peaks. Since T c is well above room 
temperature, the two-peak structure is the only one expected in the experiments of the 
following section. The breadth of the peaks is determined by the inverse coherence length 
At hp the division of spectral weight between peaks is even; for k- vectors deeper 
in the valence band more weight is found in the lower peak. This may be compared 
directly with coherence factor effects in the BCS-like phase of the model below T 3D . The 
opening of the quasigap in the density of states is visible then in A(k, u) as a splitting of 
spectral weight into two broad peaks, and the transition from metal to insulator marked 
by the progressive narrowing of these peaks until they become delta functions separated 
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by a real gap at T c . Tight bunching of peaks near the edges of the gap lead to the 
pronounced rise in the density of states there; this will again resolve into a square-root 
divergence below Tj? D . 



4. Angle Resolved Photoemission 

Quasi-one-dimensional materials were chosen for angle-resolved photoemission exper- 
iments on account of the interesting phase transitions which they undergo, and also 
because the reduced dimensional nature of their Brillouin zones simplified the interpre- 
tation of the spectra obtained. The existence of many exactly soluble models of one- 
dimensional systems makes the exploration of quasi-one-dimensional materials equally 
appealing from a theoretical point of view. The particular hope of finding evidence 
for the existence of a Luttinger Liquid has motivated a large number of photoemission 
studies of quasi-one-dimensional conductors. Among these, the most detailed studies 



have been performed on (TaSe^I, [0, f[T|, HHL an< ^ ^he discussion of the data on this 
compound is our object in this section. 

While the authors of the papers on (TaSe 2 )4/ often disagree in the detailed 
interpretation of their data, the lack of spectral weight at the Fermi energy is a 
universally observed trend. As this is one of the signal features of a Luttinger 
liquid, several researchers have concluded that (TaSe^/ and its sister compounds are 
Luttinger liquids in their conducting phase. The nature of the loss of spectral weight 
may be best probed by angle resolved photoemission, since this offers some hope of 
establishing whither the missing weight has moved. 

The quantity measured by an ideal angle-resolved photoemission experiment at 
zero temperature is the ground state electronic spectral function of the material being 
probed. This is formally equivalent to the imaginary part of the system's retarded 
Green's function, and so confers complete knowledge of its single particle properties. 
Real photoemission experiments suffer limitations of finite temperature and resolution, 
and can only probe the properties of the surface layer of material from which the 
electron emerges. Spectra also contain a substantial systematic background, which is not 
generally well understood, and which must be subtracted according to some prescription 
before the data can be fully analysed. 

Recent angle-resolved studies of the conduction band in (TaS^)^) have been 
performed above and below its CDW transition temperature, and the spectra possess 
a number of interesting features. The one most stressed in the literature is the lack 
of spectral weight at and near the Fermi energy. Related to this is the lack of a 
Fermi step in the background of spectra takena&ove the transition temperature. These 
features represent a marked departure from the usual character of metallic conduction, 
as observed in similar experiments on three-dimensional systems. We shall see below 
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that there are several other characteristics of the experimental results that are similarly 
anomalous. 

The simplest result (non-interacting electrons) for energy-distribution curves 
(EDC's) in ARPES done for initial electron momentum near the Fermi energy would be 
a delta-function. Finite experimental resolution will of course broaden this peak. We 
have convoluted our results with a Gaussian resolution function. The width is equal 
to the published estimated resolution of 160 meV. One source of (not very interesting) 
background is that of secondary electrons. The intensity of this rises steeply as the 
detected energy decreases, starting a few eV below the Fermi energy. It was found 
that, for each spectrum, a Gaussian with width of similar order to that of the band and 
centred at an energy below the band minimum could be chosen to closely mimic this 
contribution, and was duly subtracted. The overlap of the Gaussian with the interesting 
structures near the Fermi energy is small, but due to its width not quite negligible. 

As a rule, however, the most striking difference between observed EDC's and the 
ideal result is the existence of a structure resembling a Fermi edge in addition to the 
expected peak. This should be attributed to the existence of quasi-elastic scattering 
of electrons on exit from the sample, probably arising from disorder near the surface. 
It is important to remember that electrons detected in ARPES originate from near 
the surface and few surfaces are atomically flat. To take this into account, it is 
most reasonable to suppose that some fraction of the electrons have their momenta 
randomized on exit. 

The resulting spectrum is a linear combination of a true angle-resolved spectrum, 
(the spectral function), and an angle- integrated spectrum (the density of states). In 
theoretical terms, this means a combination of the imaginary part of the one-particle 
Greeen's function and the imaginary part of its trace. This picture is confirmed by the 
observation that ARPES EDC's and angle-integrated spectra taken on a single sample 
do differ only by the peak structure in the former. The relative weight in the two 
components must be determined by a fit. 

The experiments on {TaSe^zI fit this picture, with one exception. In contrast to 
experiments on three-dimensional metals, the angle-integrated spectrum near the Fermi 
energy does not resemble the expected Fermi function. Instead, the occupation falls off 
smoothly in the neighborhood of the chemical potential. Qualitatively, this may be the 
result of Luttinger liquid behavior. However, it may also result from the pseudogap in 
the LRA theory. Only a quantitative comparison can distinguish these alternatives. 

We show the comparison of experiment and theory in Fig. 6. Each plot is taken 
at a fixed angle of outcoming electron, and the inferred wavevectors are as shown. 
The data are taken from Ref. [^T|. The choice of wavevectors shown was dictated by 



the availability of experimental data. The background, fit by a Gaussian as mentioned 
above, has already been subtracted from the experimental curves. The theoretical curves 
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are plotted from Eq. 16, broadened by convolution with the Gaussian resolution function. 
Two free parameters were retained for the fit: the mean-field transition temperature 
with a best-fit value T c = 892A" and the Fermi velocity vf = 6.5 x 10 5 ms _1 , which is 
fixed by the overall linear dispersion. The parameter kp = 0.27A -1 is fixed by the fact 
that the band is one quarter filled. This gives a zero temperature gap Ao = 0.52eV 
by extrapolation of the dispersion to low temperatures. T c compares well with the 
expectation that, at a temperature around T c /4, the actual transition should take place 
and the resistivity should become activated. In (TaS , e 4 ) 2 / this occurs at T C 3D = 263K. 

Each plot shows a peak broader than the experimental resolution of 160 meV. 
The fits have been made with the stated resolution. It is evident that extremely 
good fits could be made either by increasing the width of the convolving function by 
about 50 meV (30%), by assuming some scattering of the electrons as they exit the 
material, or by assuming that the material is impure to begin with, which would give 
a momentum-independent additional width to the spectral function. Since all three of 
these alternatives involve the introduction of an ad hoc parameter, we have preferred 
to leave the theoretical curves as shown and merely note that it would be somewhat 
surprising if there was no source of broadening beyond the instrumental resolution and 
that in Eq. 16. 

Evaluating the comparisons in Fig. 6, we may say that the peak positions are given 
very well. The worst case is Fig. 6c, where the theoretical prediction is too low by 
perhaps 30 meV, and the other discrepancies are smaller. The widths are too large by 
about 20% in all cases, suggesting some relatively minor additional systematic effect. 

The momentum- integrated spectrum in Fig. 7 is obtained by integrating Eq. 16 
over k. It may be compared with the results of a angle- integrated experiments [JTTJ, 
and is clearly in good qualitative agreement with these, although the limited validity of 
our linear approximation to the free electron dispersion renders quantative comparision 
away from the Fermi energy impossible. Most significant is the movement of weight away 
from the Fermi energy, signature of the pseudo gap caused by the charge fluctuations. 

5. Fermi liquid, Luttinger liquid, or LRA liquid ? 

Four features of the LRA theory and the data are striking. 

(1) The movement of the peak position as a function of momentum is very small near 
hp. In fact, the dispersion relation, if it is defined as the peak position, apparently nears 
a quadratic maximum at k p . 

(2) The peaks are broad and symmetric. The widths are not very momentum-dependent, 
ranging only from about 400-600 meV in the range under study. 

(3) There is a strong pseudogap at all momenta, with the weight of the spectral function 
at the Fermi energy small. In addition, there is clear evidence of an energy scale 
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associated with the gap structure. This is best read off from the peak position in 
Fig. 7 as about 500 meV. This is related to the zero-temperature CDW gap in the LRA 
theory. 

(4) There are what may be called "shadow bands" . These are electronic states, or rather 
peaks in the spectral function, where no bands should be in a free-electron picture. In 
the data, a clear peak is seen even at momenta \k\ > kp. These peaks shadow the 
ordinary band in the range \k\ < kp - they are translates of the ordinary peaks through 
±2k F . 

On all four points the data agree qualitatively, and even semiquantitatively, with 
the LRA theory. 

Let us first compare these findings with the expectations of Fermi liquid theory. The 
spectral function should approach 

A(k,u) ~ 5(uj - v F (k - k F )) (18) 

as k approaches kp. Further from k F , we expect broadening due to interactions 
proportional to u 2 . Specifically: 

(1) ' In a Fermi liquid, the peak should disperse linearly through the Fermi momentum. 

(2) ' The peak should be symmetric, and the width should be resolution-limited at kp 
and broaden away from kp. 

(3) ' There is no gap or pseudogap and the only energy scales are ep > leV and 
k B T « 30meV. 

(4) ' There are no peaks when k > kp. 

It is clear, comparing points (1)' to (4)' to experiment (1) to (4), that this simple Fermi 
liquid behavior is not at all consistent with the observations. 

In the Luttinger liquid, we have a very different form for the spectral function at low 
energies. In the spinless case |[23|| , the delta function characteristic of the Fermi liquid 
is replaced by a power law singularity: 

A(k,u) ~ e{u + dF\{k-kF)\){-U + VF{k-kp)y- 1 {-UJ-dF{k-kF)y. {19) 

In this formula, denotes the step function: 0(x) =0 for x < 0, and 0(x) = 1 for 
x > 0. Two is the energy measured relative to the chemical potential so that u < 
for initial electron energies less than the chemical potential. 7 is the coupling strength 
for the electron-electron interaction. For short (finite) range interactions we expect 
< 7 < 1, (the Infinite U Hubbard model has 7 = | pi[| ). Here vp is the excitation 
velocity, which includes contributions from the kinetic energy and the interaction. The 
integrated spectral function (density of states) has the low-energy form: 

/ A{k,u)dk ~ |w| 27 . (20) 

In the case of electrons with spin, there are generally two singularities one associated 
with the charge excitations at u — VF fi {k — kp) and one associated with the spin 
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excitations at uj = vp tS (k — kp). If the velocities of the two sorts of excitations are 
very similar, then at finite resolution it may be difficult to distinguish this case from 
the spinless case. We may now compare the Luttinger liquid scenario to experiment. 

(1) " In a Luttinger liquid, for 7 < 1 the singularity or singularities disperse linearly 
through the Fermi energy. For 7 > 1, the dispersing structure becomes very diffuse near 
the Fermi energy. 

(2) " The widths that should be observed in a real experiment are not universal - they 
depend on the details of the interaction [^J. One would expect some asymmetry in the 
peaks, however. 

(3) " There is a pseudogap-like feature in the density of states. However, the exponent 
needed to fit experiment is larger than expected There is no obvious energy scale 
in the theory for repulsive interactions, though a peak may be observed in the density 
of states for |27j attractive interactions. 

(4) " There can be peaks for uj < and k > k F . However, they would be expected to be 
rather insignificant if < 7 < 1. 

Here it appears that (1)", and probably the asymmetry in (2)" and the exponent in 
(3)" are serious problems for the theory. 

We conclude that the ARPES data on (TaSe^/ are not consistent with either a 
Fermi liquid picture or a Luttinger liquid picture. 

The most decisive observation is the failure of the peaks in the EDO's to disperse 
through the Fermi energy. The shapes and the widths of the peaks as a function 
of momentum are also difficult to reconcile with these theories. The picture of very 
strong charge density fluctuations at 2k p appears, in contrast, to offer a consistent 
interpretation of all the data. We note, however, that the value of the CDW gap extracted 
from photoemission is larger than that obatined from optical conductivity or resistivity 
measurements. 

Efforts have been made in the last few years to obtain spectral predictions for the 



Luther-Emery model ||28|| . This provides a more natural framework for consideration 
of the materials in question, since CDW fluctuations are not expected to coexist with 
a Luttinger Liquid. There is not yet, however, a universally agreed prediction which 
may be compared with experiment, and the lack of a simple energy scale in the model 
makes it difficult to extract quantitaive information from the real spectra. It has also 
been brought to our attention since completion of this work, that a more complete and 
self-consistent treatment of electron-phonon interaction has recently been developed by 
McKenzie p9| . A striking feature of his theory is the asymmetry of the dispersing peaks 
in the spectral function, which is almost redolent of the Luttinger Liquid, although 
quite different in its physical origin. We have not attempted to make a quantitative 
comparison of McKenzie's spectral function with the photoemission data for (TaiSe^/, 
but note that there is no observable asymmetry in the data. 



15 



It must be stressed that these conclusions apply to this particular system only; 
experiments performed on other quasi-one-dimensional systems suggest that the relative 
importance of electron-electron and electron-phonon interaction effects is strongly 
dependant on the system in question. Even for (TaSe^/, correlation effects may 
be masked by the strong charge fluctuations only within a particular parameter regime; 
it is for example possible that the application of pressure to the chains would relatively 
strengthen the interaction effects. In other materials, and particularly in some of the 
quasi one- dimensional SDW compounds, there would appear to be evidence for strongly 
one-dimensional correlation effects. 
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Figure captions 



Figure 1. Phases of a Peierls insulator - from Kohn Metal with phonon frequency 
Q 2 kf to Band Insulator with gap A(T). 



Figure 2. Schematic structure of the Ta chain in (TaSe^-f • 



Figure 3. Feynman diagram for the electron self-energy in the LRA theory, with the 
phonon self-energy indicated by the box. 



Figure 4. LRA density of states as a function of energy and temperature for the 
conducting phase. 
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Figure 5. Plot of n(k), the occupation number, for temperatures above, below, and 
at (solid line) the 3D ordering temperature. 



Figure 6. Single parameter fits to ARPES spectra taken on (TaSe^-f in conducting 
phase, for various wavevectors for which data were available. 



Figure 7. Momentum integrated spectrum found from LRA model in conducting 
phase, showing suppression of the density of states at the chemical potential /i = 0. 




Kohn Metal 
^2k f (T) 



d z orbital 



Ta 

Se 



, Se Rectangles 

"O are rotated with 
regard to one 
^ ^ another. 



A 



Direction of Ta chain 



a) k = k f 




-1.2 -1 -0.8 -0.6 -0.4 -0.2 

Binding Energy in eV 



b) k= 1.19 k f 





Binding Energy in eV 



d ) k = 1 .03 kf 




-5 l= 

-1.2 -1 -0.8 -0.6 -0.4 -0.2 



Binding Energy in eV 



e) k = -0.89k f 



Experiment o 
/ \ Theory 



y° o o \ 



O v / 
/ 

O / 
O I 

I 

' / 



° / 

/ 

/ 

/ 

O O / 

o 



\ 



\0 
\ 

\ o 
\ o 



\ 



\ o 



o 

\ o 



^ o * 
"V o 



-1.2 



-0.8 



-0.6 



-0.4 



-0.2 







Binding Energy in eV 




Binding Energy in eV 




Energy in eV 



